Macroscopic Resonant Tunneling in the Presence of Low Frequency Noise 
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We develop a theory of macroscopic resonant tunneling of flux in a double-well potential in the 
presence of realistic flux noise with significant low- frequency component. The rate of incoherent flux 
tunneling between the wells exhibits resonant peaks, the shape and position of which reflect qualita- 
tive features of the noise, and can thus serve as a diagnostic tool for studying the low-frequency flux 
noise in SQUID qubits. We show, in particular, that the noise-induced renormalization of the first 
resonant peak provides direct information on the temperature of the noise source and the strength 
of its quantum component. 



Superconducting qubits based on different forms of 
quantum dynamics of magnetic flux in SQUIDs continue 
to demonstrate steady progress - see, e.g., [H, 0, 01- Nev- 
ertheless, the low-frequency flux noise in SQUID struc- 
tures still provides the major obstacle to their further de- 
velopment to the level necessary for quantum computing 
applications. Because of this noise, quantum coherence 
is essentially reduced to the states with the same aver- 
age flux (e.g., to qubit operation at the "optimal point") 
Only very limited coherence exists between the 
states which differ by their flux values In supercon- 
ducting charge qubits, the low- frequency charge noise is 
most probably produced by elementary charges tunneling 
between impuritj^states localized in the insulator parts 
of the structure [5|. In contrast to this, it is more diffi- 
cult to construct a possible model of the low-frequency 
noise of magnetic flux. Elementary magnetic moments 
of nuclei or electrons [6| are small and require very large 
concentration of impurity states to explain the typical 
value 6^ ~ 0.1 m"I>o of the observed noise 0, So the 
origin of the low-frequency flux noise in SQUID struc- 
tures remains not fully understood, making it interesting 
to probe its properties in different experiments. In this 
work, we show that one of the interesting dynamic pro- 
cesses in SQUIDs, macroscopic resonant tunnelin g of flux 
between the wells of the double-well potential [ol Il0l| can 
serve as noise diagnostic tool. In particular, the shape 
and position of the resonant flux tunneling peaks reflects 
the temperature and the strength of quantum component 
of the low-frequency noise. 

Specifically, we consider a SQUID in the regime when 
its flux potential has the double- well form (Fig. [T|) . The 
extraction of precise parameters of the system, e.g., the 
spectrum ej of the energy levels in the two wells, re- 
quires its direct numerical simulation. Considerable in- 
sight into the flux dynamics in this system can still be 
gained from analytical treatment based on the general 
features of the spectrum. The states localized in each 
well are separated by large energy intervals on the or- 
der of the plasma frequency ujp, while the tunnel cou- 
pling A between the states in the opposite wells is small 
(?i = fcs = 1), A ^ Wp. This means that the rate of flux 



FIG. 1: Schematic energy diagram of the macroscopic res- 
onant tunneling in the double-well system between (a) two 
lowest energy states and (b) two excited states in each well. 



transfer between the two wells exhibits resonant peaks 
as a function of bias between the wells [§| whenever the 
energy distance e between a pair of levels in the opposite 
wells becomes small, e <^ ojp. In the vicinity of such a 
resonance, the influence of the environment on the flux 
dynamics in this system of levels separates into two dif- 
ferent effects [l3| even if the dissipation is produced by 
physically one and the same environment. The flrst is 
"intrawell" relaxation that causes the dissipative transi- 
tions within each well, and is sensitive to the properties 
of environment at large energies LOp, and the second is 
"interwell" relaxation associated with relaxation dynam- 
ics within the two tunnel-coupled states, and sensitive 
to the environment properties at low energies ^ LUp. In 
terms of the coupling of flux to the environment, inter- 
well relaxation is dominantly affected by fluctuations of 
the bias e that shift one well of the double-well system 
as a whole relative to the other. The purpose of this 
work is to extend the previous theories [lO|, 



III of the 



resonant flux tunneling to the regime of strong interwell 
relaxation. This regime is made experimentally relevant 
by apparently unavoidable low-frequency flux noise in the 
SQUID structures. The developed approach should also 
be valid for the description of tunneling in other double- 
well systems besides flux qubits. 

We start by calculating the transition rate for tunnel- 
ing between the lowest energy levels |0) and |1) in the 
"left" and "right" wells, respectively (Fig. [1^). Such a 
process can be described by a two-state model: H = 
Hs + Hb + -ffint, where Hb is the environment's Hamil- 
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tonian and 



^^int — —-^O-zQ 



(1) 



are the system and interaction Hamiltonians respectively. 
Here, ax,z are the Pauh matrices and Q is an operator 
acting on the environment. We adopt the Gaussian ap- 
proximation for which the dissipative dynamics is char- 
acterized completely by the spectral properties of the 
reservoir force Q, hence the reservoir Hamiltonian Hb 
does not need to be made explicit. If the environment 
is in equilibrium at temperature T, the (unsymmetrized) 
spectral density of the noise Q, given by 



S{co) 



1 

2^ 



dt e'^'{Q{t)Q{0)) 



where (...) = Trsjp...}, can be expressed in the usual 
way in the basis of energy eigenstates |n) of Hb with 
eigenvalues £„ and equilibrium probability density pn 



=^Pn\{n\Q\m)\'^S{£n - £„ 



(2) 



of in general strong interaction with the low-frequency 
environment, we extended the starting time of the evo- 
lution U{t) to —oo. Substituting ([4]) into ([3]), we find 

ro^i(e)-ri^o(-e) = r(e), with 

A 2 roc 

r(e) = — Rey^ dte-*([/l(i)t/+(t)C/t(0)C/_(0)), (5) 

where U± — Texp{±(i/2) f^^Q{T)dT}. Here, we have 
summed the rates ^i^f over the initial (with equilibrium 
density matrix p) and final states of the environment. 
The correlator in Eq. ([5]) can be calculated in the Gaus- 
sian approximation by expanding each of the operators 
U± up to the second order in Q, averaging the result and 
exponentiating it back. This gives: 



exp 



t M 
dT / dT'(Q{T)Q{T')) 
J -oo 



Expressing the noise correlator in this equation through 
the spectral density ([2]), we reduce Eq. ^ to 



We find the rate T of flux tunneling between the reso- 
nant states |0) and |1) in the incoherent regime A ^ W, 
where W is the noise-induced resonance width. In this 
case, r can be calculated in the lowest-order perturbation 
theory in the tunneling term V = Aax:/2 in the Hamil- 
tonian IT]) . We assume that the width W is still small on 
the scale of the plasma energy cOp, so that the system dy- 
namics can be reduced to the two resonant levels even in 
the presence of noise. This condition is satisfied in typi- 
cal experiments, where W is on the order of a few GHz, 
while LOp ~ 30 GHz [^, . The lowest-order transition 
rate between the states \i) and |/) can be obtained using 
Fermi Golden rule: 

T,^f = 2n\mV\f)\'S{E,-Ef) 



dte 



UEi~Ef)t 



\mf)\ 



= 2Rc / dt{i,t\V\f,t){f,{)\V\i,Q). (3) 
Jo 

where Eij are the eigenvalues of the unperturbed total 
(system -I- environment) Hamiltonian for states |j),|/). 
In our case, these states are the systems resonant flux 
states |0) and |1) plus the states of the environment at 
the beginning and the end of the transition. The time 
evolution of these states can be written as {a — i, f): 



|a,t) = t/(t)e*^"^*/2|a,0). 



(4) 



Here, U{t) = T exp{(i/2) cr^ /^^ Q(r)(ir} is the evolu- 
tion operator in the interaction representation, with T 
denoting the time ordering and Q{t) = e^^^'Qe"*^^*. 
To avoid description of the process of "switching on" 
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— Re 
2 JO 



dte'"' exp <^ / dujS{uj 



(6) 

The integral over time in Eq. ^ converges on the scale 
which defines the width W of the resonance. The 
line-shape of the resonance is determined therefore by the 
noise frequencies uj < W. If S'(w) is essentially constant 
in this frequency range, S{uj) — S{0), i.e., if the noise is 
sufficiently broad-band and/or weak, the tunneling rate 
([6]) has Lorentzian line-shape: 



2 e2 + 1^2 ' 



W = 7rS'(0) 



(7) 



In general, Eq. © can lead to line-shapes that are 
not Lorentzian. If the cut-off frequency Uc of the low- 
frequency part of the noise satisfies the condition ujc ^ 
W (i.e., the noise is strong and has low- frequency) , then 
one can expand the exponent e~*"* up to the second order 
in cot to obtain: 



2W^ 



= / dujSico), 



ep^V 



duj 



Siu;) 



(8) 
(9) 



We assume that both integrals in ([9]) are finite, with the 
second integral being understood as a principle value. 
We see that, as usual, the noise with a significant low- 
frequency part yields a Gaussian line-shape, in contrast 
to the Lorentzian line-shape ([7]) in the case of the broad- 
band noise. Equation ([8]) shows also that, in general, 
the low-frequency noise not only broadens the resonance 
but also shifts it to the non- vanishing bias e ~ e^. This 
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means that effectively the first resonant peak spUts in 
two, described by r(e) and r(— e), for the two different 
directions of tunneling. 

If the environmental source of the low-frequency noise 
is in equilibrium at temperature T, then symmetric and 
antisymmetric (in frequency) parts of the noise intensity, 
S±{lj)=[S{uj)±S{—lj)]/2, are related by the fluctuation- 
dissipation theorem: S-{uj) — iS'+(ti;) tanh(w/2T). The 
antisymmetric noise gives also the dissipative part 
a{LO,T) of the hnear response (e.g., conductance) of the 
environment generating the noise: S-{uj) = uja{uj,T). 
Since in ([9]) the width W and shift Cp of the Gaussian tun- 
neling peak are determined, respectively, by symmetric 
and antisymmetric parts of the noise: = J dujS^{LL!) 
and Ep — J duj S - (oj) / u! , the fluctuation-dissipation the- 
orem relates W and ep. This relation is particularly 
simple in the case of low-frequency noise, when all the 
relevant frequencies, uj < tOc, are small on the scale of 
temperature T: ui <^ T. For typical cut-off frequen- 
cies on the order GHz, this condition holds even at 
sub-Kelvin temperatures characteristic for experiments 
with solid-state qubits (1 GHz ~ 50 mK). In this case, 
tanh(tj/2T) ~ uj/2T, yielding 

= 2Tep . (10) 

Equation (fTO)l is one of the main conclusions of the the- 
ory developed in this work. It shows that, qualitatively, 
non-vanishing splitting 2ep of the two parts of the first 
resonant peak is a manifestation of the finite quantum 
component of the low- frequency flux noise. In the regime 
of classical noise (which is valid, due to the fluctuation- 
dissipation theorem, only as a high-temperature approxi- 
mation), the spectrum S{uj) is symmetric in w and ep — 0. 
Quantitatively, this classical limit is reached at T ^ W , 
when ep<SLW and the two peaks merge. 

The above method of the calculation of the tunneling 
rates can be generalized to the case of resonant tunnel- 
ing between two higher energy levels in the wells. In that 
case, in addition to the interwell relaxation process, we 
will also have intrawell relaxation. We still use |0) and |1) 
to denote the two resonant flux states, which can now be 
excited states in their corresponding wells (see Fig. [T|d). 
In addition to Hamiltonians ([1]) , we need to add to the to- 
tal Hamiltonian two other terms: Hq = X^fc^^^o i 
that describes the energy of the flux energy eigenstates 
|fc) with /s 7^ 0, 1, and 

llr= ^ <^fefeK|fc)(fc'l + |fc')(fc|)Q- (11) 

that produces the intrawell relaxation. The sum in 
goes over the pairs of states and that belong to 
the same well and (/>fcfc' are the coupling matrix elements 
between these states that drive the relaxation transitions. 

We make two main assumptions about the properties 
of decoherence. First is that the high-frequency compo- 
nent of the noise is sufficiently weak, so that the intrawell 



relaxation rates in the two wells are small on the scale 
of Wp, Vh,R ^ f^p, so that it makes sense to discuss res- 
onances in the tunneling rate. The relaxation rates can 
still be large on the scale of tunneling strength A. Un- 
der a natural assumption T <C Wp, effects of the weak 
(on the scale of Wp) intrawell relaxation depend only on 
the spectral components of the w ~ Wp that enter into 
the relaxation rates Fi^^j, and we can characterize them 
directly by these rates. 

We find the rate F of flux tunneling between the res- 
onant states |0) and |1) the same way as before. Since 
the environmental modes that contribute to the inter- 
well and intrawell relaxations are from different parts of 
the environment spectrum and the coupling Hamiltoni- 
ans iJint and Hr commute (the matrix elements <\)kki are 
non-vanishing only between the states within the same 
well), Eq. ([4]) can be written as direct combination of 
the evolution due to the high-frequency noise U(t) = 
Texp{— i Hr{T)dT} and the low-frequency noise U{t), 
as defined below (H]): 

\a,t) = U{t)U{t)e""''*\a,0), 

where Hrit) = e'^"^+"'>'>* Hre-'^"^+""'>K Since the re- 
laxation is assumed weak on the scale of the plasma fre- 
quency, such a trace over the high-frequency part of the 
environment can be done in the relaxation approxima- 
tion, when the relaxation dynamics is determined by the 
transition of the lowest (second) order in that are 
characterized by the relaxation rates T^n. Then, 

TTB{p{j\UHm{j'\Um')}-e-''\ (12) 

where j, j' denote the two opposite resonant flux states 
coupled by tmmcling: j, j' G {0,1}, j ^ j' , and 7 = 
(Fi-|-F/j)/2. Equation (fT2|) is written in the form conve- 
nient for the calculation of the tunneling rate in Eq. 
It can be understood more directly if one notices that 
it gives the time evolution of the off-diagonal element 
of 11^6 flux density matrix averaged over the en- 
vironment. According to the standard theory of weak 
relaxation (see, e.g., [13]), such a matrix element decays 
with the rate equal to the half-sum of all the relaxation 
rates out of the states \j) and |/) - cf. Eq. (|12p . After 
the high- frequency part of the environment is traced out, 
expression for the flux tunneling rate takes the form 

F(e) = —Re / die^*^"^)* exp <^ / dujS{uj)—-^ — \ . 

(13) 

This equation describes higher resonant peaks, when 
the flux tunnels from the ground state in the initial (e.g., 
left) well into an excited state of the target (right) well. 
In this case, 7 — Vr/2, and the line-shape of the reso- 
nance ()13|) is determined by the combined action of re- 
laxation broadening and the low-frequency noise. If the 
noise is strong, using the same approximation as for the 
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first resonance, we get: 



«7 



(14) 



where the parameters are given as before by Eq. ^ , and 
w{x) is the complex error function defined as 



w(x) = e [1 - erff-ij:)] = e "^'^ 



Equation (fT4|) describes the transition between the Gaus- 
sian and Lorentzian hne-shapes in the noise- (W ^ 7) 
and the relaxation-dominated (W ^ 7) regimes, respec- 
tively. Independently of this, if the noise is quantum and 

7^ 0, the higher resonant peaks are also shifted by 
away from their positions determined by the energy lev- 
els in the well that follow from the Schrodinger equation 
for the flux dynamics in the absence of environment. 

Finally, we consider the microwave-activated tunneling 
between the two wells. For the simplest case where the 
tunneling occurs between the two lowest energy levels of 
each well, the process can be described by the Hamil- 
tonian H^f — —eaz + V cos fit, where V is an operator 
that couples microwaves to the flux dynamics. Unitary 
transformation U = e-'l^Vz)^^ changes this Hamilto- 
nian into Hrf = UHsW + iUW, and gives in the rotat- 
ing wave approximation: Hrf ~ ^\{'^<^z + V<^x), where 
e — e~D,<$^D,is the detuning of microwave radiation 
and V — (0|y|l). In the absence of the environment, the 
system would undergo Rabi oscillations with the Rabi 



frequency fifl = y -|- jVp. If coupling to the environ- 
ment is strong enough so that V ^ W, then the tunnel- 
ing becomes incoherent with the rate given by ([5]), but 
with e and A replaced by e and \V\, respectively. 

In another interesting regime, the microwave excites 
the system from the ground state |fc) to an excited state 
|0) in the first well, and it tunnels then to state |1) in the 
target well (Fig. [T]). We make use of the same assump- 
tion of small tunnel coupling A that can be treated as 
perturbation. In addition to previously used condition 
A ^ W, this assumption implies also that the tunnel- 
ing rate r(e) is small compared to the rates of excita- 
tion/relaxation processes within each well. The balance 
between excitation and relaxation (Tl) in the well es- 
tablishes the stationary occupation po of the state |0). 
The standard solution for the relaxation dynamics of 
the density matrix of the system (see, e.g., |lO|) gives 
po = |l^lV(e' + 2|^P+ri/4).^ Here, again, e~ = eo-efc-l^ 
is the detuning energy, and V — {0\V\k). At this stage, 
one can repeat the same steps in the calculation of the in- 
coherent tunneling rate between the states |0) and |1) as 
above, and obtain the rate F = poT, where F is given by 
()14|) . but now with the relaxation broadening determined 
by the relaxation in both wells, 7 = (F^ -|- Fj^)/2. Quali- 
tatively, as before, the shape of the resonant peak of F is 



determined by a convolution of the Lorentzian broaden- 
ing due to intrawell relaxation and Gaussian broadening 
due to interwell noise. 

In conclusion, we have developed the theory of macro- 
scopic resonant tunneling in flux qubits in the presence of 
Gaussian low-frequency flux noise. The tunneling rate is 
given in general by Eq. ([13]) . In the case of strong noise 
and tunneling between the two lowest energy levels in 
each well, the resonant peaks have Gaussian shape ([8]). 
If the noise source is in equilibrium at temperature T, 
there is a fundamental relation pH]) between the width 
W of the resonant peaks and shift Cp of their position 
in energy. The shift ep provides direct measure of the 
strength of the quantum component in the flux noise, 
and leads to splitting of the first resonant peak of flux 
tunneling. Some of the predictions of the present work 
has been already tested in experiment 1^, which shows, 
in particular, that Eq. (fTO|) is a convenient tool to deter- 
mine whether the low-frequency flux noise is produced 
by an equilibrium quantum source. More experiments, 
especially with microwave excitations, can be performed 
to further test our theory. 
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